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Abstract. The GiBUU model, which implements all reaction channels relevant at medium 
neutrino energy, is used to investigate the neutrino and antineutrino scattering on iron. Results 
for integrated cross sections are compared with NOMAD and MINOS data. It is shown, that 
final state interaction can noticeably change the spectra of the outgoing hadrons. Predictions 
for the Minerz/a experiment are made for pion spectra, averaged over NuMI neutrino and 
antineutrino fluxes. 
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1. Introduction 

Neutrino and antineutrino scattering on nuclei for neutrino energies above 30 GeV was studied in 
several experiments starting from the 80s. Theoretically they were successfully described within 
the quark parton model as Deep Inelastic Scattering (DIS) processes. Recent measurements 
by MINOS and NOMAD colalb orations also covered the intermediate energy region. Here the 
neutrino reactions are not so easy to model, because of the overlapping contributions from QE 
scattering, resonance production and background processes. This requires complex approaches 
that take all of the relevant channels into account. 

Nuclear effects in neutrino reactions can also be studied in detail nowadays. The Minerz/a 
experiment intends to perform measurements on Plastic (CH), Iron, Lead, Carbon, Water and 
liquid Helium targets in the NuMI beam, which would directly allow to compare nuclear effects 
on various nuclei. Besides muon detection, this experiment will also be able to resolve various 
final states by identifying the tracks of the outgoing hadrons. Theoretical modelling of such 
exclusive reactions requires realistic approaches to the desctiption of the initial and final state 
interactions in target nuclei. 

Both these requirements are satisfied by the GiBUU transport model [U [2], which we use 
here to study neutrino and antineutrino scattering on iron. Our results are compared with the 
recent MINOS and NOMAD data. Predictions are also made for the spectra of the outgoing 
particles. The calculations are done without any fine tuning to the data covered here with the 
default parameters as used in the GiBUU framework. 



2. GiBUU transport model 

The GiBUU model was initially developed as a transport model for nucleon-, nucleus-, pion-, 
and electron- induced reactions from some hundreds MeV up to tens of GeV. Several years 
ago neutrino-induced interactions were also implemented for the energies up to a few GeV. 
Recently the code was extended to describe also the DIS processes in neutrino reactions. Thus, 



we can study all kind of elementary collisions on all kind of nuclei within a unified framework. 
The model is based on well-founded theoretical ingredients and has been tested against various 
nuclear reactions. For a detailed review of the GiBUU model see [2j. 

GiBUU describes all processes relevant at medium energies, the cross section is calculated as 
^tot — '^QE + (^RES + (^BG + (^Dis- Our approach to quasielastic (QE) scattering, resonance 
(RES) production and background (BG) processes is described in [3l[l]. The DIS scattering is 
included as Pythia simulation. 

In the region of the shallow inelastic scattering, that is at moderate invariant masses, 
1.6 GeV < W < 2.0 GeV, there is a potential problem of double counting. Here the same 
physical events can be considered as originating from decays of high mass baryonic resonances 
or from DIS. In the GiBUU code we use the ansatz, that both RES, BG and DIS processes 
contribute in this region. While the RES and BG contributions are smoothly switched off in this 
region and DIS contribution is switched on. With this choice, the DIS events become noticeable 
at neutrino energies around 3 GeV. In essence, the DIS processes below at lower W account 
for the resonances whose electromagnetic properties are not known and for the non-resonant 
processes giving several mesons in the final state beyond the one pion background. 

3. Integrated cross sections. 

As was measured in 80s, the DIS cross section grows linearly with energy. So at high 
neutrino energies the data are conveniently presented as cross section per energy a^Q^/E,j. 
Despite the measurements being made on nuclear targets, the world average values are given 
^ for isoscalar target: atot/E^lu) = 0.667 ± 0.014 • 10~^^ cm^/ GeV for neutrinos and 
atot/Ejj{iy) = 0.334 lb 0.008 • 10~^^ cm^/ GeV for antineutrinos. Such an approach is meaningful, 
only if nuclear corrections are neglected. 
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Figure 1. (Color online) Total cross section per nucleon for neutrino (left) and antineutrino 
(right) induced reactions on iron target. Various prescriptions to include nuclear effects in DIS 
are compared. MINOS data (full circles) are from [SIE], IHEP-JINR data (full diamonds) are 
from [8J. 

The actual value of such corrections for neutrino reactions is not known so far, because of both 
experimental inaccuracies and difficulties in the theoretical description. On one hand, nuclear 
parton distributions, based on electromagnetic scattering data and intended for description of 
both charged lepton and neutrino reactions, were introduced. For a review and a list of recent 
parametrization see, for example, [9]. On the other hand, recent investigation [10^ 111] showed, 
that in neutrino reactions nuclear corrections to parton distributions are at the same level as for 
electrons, but have a very different dependence on the Bjorken x variable. The topic remains 



controversial, with the hope that future precise Mineri^a results on various targets will clarify 
the situation. 

As we already mentioned, the GiBUU code uses PYTHIA for the simulation of DIS processes. 
In the GiBUU simulation the neutrino interacts with one initial nucleon, bound in the hadronic 
potential and having nonzero Fermi momentum. In order to be able to use the PYTHIA 
event simulator, we have to provide some quasi-free kinematics as inputs to PYTHIA. Various 
prescriptions to do this result in a 5 — 7% difference in the results. The corresponding cross 
sections (denoted as "F-NO", "F-CM", "F-THRE") are shown in Fig [H Nuclear parton 
distribution functions from [121 ^-^so implemented as one of the options to use (to avoid 
double counting, nuclear potential and Fermi motion in such calculations are switched off). The 
result ( "nuclearPDF" ) as well as the free cross section for iron composition ("26p-|-30n") are 
also shown in Fig[TJ At the moment we consider the various prescriptions mentioned above as 
intrinsic uncertainty of the GiBUU code, reflecting the lack of our understanding the nuclear 
effects. No other event generator, as far as we know, accounts for nuclear corrections in high- 
energy neutrino reactions. 

Fig. [1] shows, that our calculations are in a good agreement with the recent neutrino data, 
which are also consistent with each other. This figure shows that the decreasing slope of our 
curves for the neutrino cross section is in agreement with that of the data. This slope is not 
taken into account in deriving the world-average value, where it was assumed to be negligible. 
For antineutrino the agreement is good for Ep > 25 GeV. For lower energies our curve is above 
the recent MINOS data, but below the IHEP-JINR results [8]. The overall agreement of our 
calculations with the data is therefore better than the agreement of the data with each other. 

4. Final state interactions and change of the final hadronic spectra. 

After being produced in the initial interaction, outgoing hadrons propagate throughout the 
nucleus. In GiBUU this process of final state interactions (FSI) is modeled by solving the semi- 
classical Boltzmann-Uehling-Uhlenbeck equation. It describes the dynamical evolution of the 
phase space density for each particle species under the influence of the mean field potential, 
introduced in the description of the initial nucleus state. Equations for various particle species 
are coupled through this mean field and also through the collision term. This term explicitly 
accounts for changes in the phase space density caused by elastic and inelastic collisions between 
particles. FSI decrease the cross sections as well as significantly modify the shapes of the final 
particle spectra. Such change was seen, for example, in photo-pion production [13j and is 
described by the GiBUU with a good accuracy. A similar change should be observed in neutrino 
reactions. 

Fig.[2]shows the vr^ (left panels), vr" (middle panels) and ir^' (right panels) spectra for neutrino 
(upper panels) and antineutrino (lower panels) NuMI fluxes for 1-pion events. The green dashed 
lines show the kinetic energy (T^r) distributions without FSI, i.e. of pions produced in the initial 
neutrino vertex. The solid red lines show the distributions after FSI, i.e. of pions that made it 
out of the nucleus. Such spectra can also be calculated for any other outgoing particles (protons, 
kaons, eta) for any predefined final state and should be measurable in Mineri/a experiment. 

For dominant channels (vr"*" production for neutrino reactions and vr" in antineutrino ones), 
the FSI decrease the cross section at T^r > 0.2 GeV. This is mainly explained by pion absorption 
through vrA^ — )• A following by AA^ — t- NN. Pion elastic scattering in the FSI also decreases the 
pion energy, thus depleting the spectra at higher energies and accumulating strength at lower 
energies. Thus, an increase of the cross sections is observed at T^r < 0.15 GeV; altogether this 
leads to a significant change of the shape of the spectra. 

Scattering can also lead to pion charge exchange. For neutrino-induced reactions, the 
TT^n ir^p scattering in the FSI is the main source of side-feeding for the vr'^ channel, leading 
to a noticeable increase of the tt^ cross section at low T^^. The inverse feeding is suppressed. 




Figure 2. (Color online) Pion kinetic energy distributions for neutrino and antineutrino induced 
reactions for 1-pion production (one pion of a given charge and no other pions are produced). 
Calculations are for the NuMI low-energy mode neutrino/antineutrino fluxes. 



because less vr than vr^ are produced at the initial vertex. The same mechanism of side feeding 
from dominant to sub-dominant channel through 7r~p — t- ir^n is working for antineutrino induced 
reactions. 

For the least dominant channel (vr^ production for neutrino reactions and vr^ in antineutrino 
ones), the FSI (in particularly side feeding) represent the main source of the events observed; 
thus a dramatic FSI effect. 
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